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Abstract 

A recent Particle Image Velocimetry (PIV) experiment in He II counterflow around a cylindrical 
obstacle showed the existence of apparently stationary normal fluid eddies both downstream (at 
the rear) and upstream (in front) of the cylinder. This rather surprising result does not have an 
analogue in experimental observations of classical fluid flows. We suggest that the explanation for 
the apparent stability of such eddies can be provided entirely from the viewpoint of classical fluid 
dynamics. We also discuss a possible connection between the emergence of the normal fluid eddies 
and the polarization of the vortex tangle in superfluid. 

PACS numbers: 67.40.Vs, 47.37.+q, 47.27.-i 
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I. MOTIVATION 



Particle Image Velocimetry (PIV), which has been, for several decades, a standard tech- 
nique of flow visualization in classical fluid dynamics, is based on tracking the motion of 
solid particles injected into the fluid. Provided the particles are sufficiently small, they can 
be expected to follow faithfully the fluid flow. 

Recently implemented in helium II ll| , the PIV technique has already yielded some non- 
trivial and unexpected result s_^ [j, 0, Q followed by the attempts of their theoretical 



interpretation (see e.g. Refs. [6|, |7 |). Among surprising experimental results is the recent 
observation [3j of the apparently stationary normal fluid eddies in the thermal counterflow 
past a cylinder. In the cited work, Zhang and Van Sciver visualized the motion of small 
particles in the thermal counterflow around the cylinder of diameter D = 0.635 cm fixed 
in the center of rectangular channel of a cross-section 3.89 x 1.95 cm 2 . The counterflow 
was produced, in two separate experiments, by the heat flux q = 0.4 and 1.12 W/cm 2 at 
temperatures T = 1.6 and 2.03 K, respectively (corresponding to the Reynolds numbers 
Re_D = pDv n /p Q = 4.1 x 10 4 and 2.1 x 10 4 , where p = 0.145 g/cm 3 is the density of the 
liquid helium, and v Q and p n are the mean velocity and the viscosity of the normal fluid). 
Solid particles used for visualization in the PIV experiments were polymer microspheres of 
diameter 1.7 /im and density 1.1 g/cm 3 . 

In these experiments Zhang and Van Sciver observed the formation of large-scale eddies of 
the particulate motion located both downstream and, surprisingly, upstream of the cylinder 
with respect to the normal flow (it should be noted that these structures were somewhat 
more discernible in the second of two experiments, i.e for Re£> = 2.1 x 10 4 , in which case 
the corresponding mean normal fluid velocity was v n ~ 2.2cm/s). These, apparently stable 
vortices of the particulate flow field were located at distances about 3 cylinder radii from 
its center at the angles ±45° and ±135° to the axis along the undisturbed flow through the 
center of the cylinder. 

In order to interpret these observations, the following question should first be addressed: 
what do tracer particles actually trace? It seems natural to assume that in the cited exper- 
iment, because the Stokes drag of small particles is much larger than other forces exerted 

n 

by the normal fluid (see e.g. Ref. [6|]), the particle trace the normal fluid. However, solid 
particles also interact strongly with quantized vortices which may reconnect to the particle 
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surface; such interactions may lead to the appearance of the additional force exerted on 
particles by the superfluid. This additional force was actually measured in another Zhang 
and Van Sciver's experiment [2] on particle sedimentation in the thermal counterflow, and 
later analyzed theoretically [7| by the authors of the present work. 

A second question is: do the circulation cells of the particulate motion map, in the 
experiment 3|, the normal eddies, or result from complex interactions of particles with both 
the normal fluid and quantized vortices in the superfluid component (in which case the 
eddies in the normal fluid may not even exist)? Since in this experiments the vortex tangle 
was relatively dilute (see more detaled discussion below in Sec. N} . we can expect with 
some confidence that the particle motion indeed maps, at least qualitatively, the normal 
flow. The observed normal flow patterns do not have a classical analogue: in classical 
fluid dynamics, stationary eddies upstream of the cylinder have never been observed, and 
downstream of the cylinder one would rather expect the Von Karman vortex street than a 
pair of apparently stationary eddies. In Ref. [3] it was suggested that the observed large 
normal vortex structures were caused by the complex interaction between the two fluid 
components of He II. 

In the present work we argue that, while the emergence of the large-scale vortex structures 
is, most likely, caused by the mutual friction between quantized vortices and the normal fluid, 
the apparent stability of the observed eddies, both at the rear and in front of the cylinder, 
can be explained entirely from the viewpoint of classical fluid dynamics without invoking an 
interaction between the normal fluid and superfluid vortices. 

Because the Reynolds numbers in the cited experiment were large (~ 10 4 ), instead of 
considering the turbulent viscous flow past a cylinder we will analyze a much simpler model 
of the plane motion of two point vortices of opposite polarities in the inviscid, potential 
flow around a disk. We shall show that there exist stationary configurations of the vortex- 
antivortex pair, both behind and in front of the disk. Such configurations are unstable, i.e. 
any perturbation of the stationary configuration leads to eventual sweeping of the vortex 
points away from the disk. However, we shall show that there exist stationary configurations 
such that, on the time scale corresponding to the duration of experiment 3J], the vortices 
located sufficiently close to the corresponding stationary points will remain in the close 
proximity of their initial locations. Moreover, such stationary points are positioned as the 
apparently stable eddies seen by Zhang and Van Sciver. In conclusion we shall also discuss 
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a possible connection between the emergence of the normal fluid eddies and the polarization 
of the vortex tangle in the superfluid component of He II. 



II. LAGRANGIAN EQUATIONS OF MOTION OF POINT VORTICES IN THE 
INVISCID FLOW AROUND THE DISK 



We consider the two-dimensional inviscid, potential flow, with velocity U at infinity, 
around a circular disk of radius a. We are concerned with a motion in such a flow of a 
vortex-antivortex pair, i.e. of two vortices of opposite polarities but the same circulation, 
T. Let the positions of the vortex points on the complex plane be Z\(t) = x\(t) +iyi(t) and 
z 2 {t) = x 2 (t) + iy 2 (t), see Fig. [TJ The reason why the vortex z\ has negative (clockwise) 
circulation, T < 0, will become clear in Sec. IIHI 

Each of these vortices moves as a fluid point in the superposition of the imposed flow 
around the disk and the flow field of another vortex. Complex potential of the superposition 
of the potential, uniform at infinity, flow around the disk and the flow created by the vortex 
point at, say, z = Z\ is js] 



w(z) = U[ z + 



+ 



ln(z — z{) — In 



a — zz\ 



: / 2iri 

where z* = x — iy denotes the complex conjugate of z = x + iy. 
It is convenient to introduce the non dimensional variables 
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so that the non-dimensional complex potential is (from now on the primes are omitted) 
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The x- and ^-components of the fluid velocity are, respectively, u = Re{dw/dz} and v = 
— Im {dw/dz}. 

The flow potential yields the following Lagrangian equations of motion of two vortex 
points, Zi(t) and -2 2 (t): 



^ = Im{/j(«i, z 2 )}, 



where j = 1, 2, and 
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(5) 



/lfe ,^ ) = 1 _^ + 2U {_i_ + _ J _i_}. (6) 

III. STATIONARY POSITIONS OF POINT VORTICES 

In relation with the experimental observations described in Ref . , we will now determine 
whether there exist stationary locations of the vortex-antivortex pair, both downstream and 
upstream in the reference frame of the disk. On the complex plane, such stationary positions 
can only be the complex conjugate points, say, z\ in the upper, and z\ = (z^)* in the lower 
half-plane, respectively. Since the imposed flow is in the positive x-direction, the velocities 
of the both vortex points can be zero only in the case where the vortices located at z\ and 
z\ have the negative and the positive polarity, respectively. This implies T < 0, so that 
the non-dimensional parameter A defined by the last relation (j2j) is positive. To ensure that 
the point vortices located at z\ and z\ are stationary, it is sufficient to require that the 
complex velocity at z% = (z®)* is zero, i.e. = fi{Z\-, (2?)*) = 0, where the function / 2 

is defined by relation (jH]). Introducing the function F(x, y) = yz*(l — (z*) 2 )f(z), we obtain 
the following pair of coupled equations for the coordinates of the stationary point: 

ImF«y°) = 0, ReF« y°) = 0, (7) 

where 

lmF{x, y) = Axy{y-\){x 2 -y 2 -l) ) (8) 

ReF(x, y) = (x 2 - y 2 )(l - x 2 + y 2 )(y - A) + Ax 2 y 2 (y - A) - y(l - x 2 + y 2 ) + 2\y 2 . (9) 

Consider first Eq. ([8]). Obviously, y° 7^ 0. For the root y® = A, from Eq. dHJ) we find 
x\ = ±yl + A 2 , so that (x°) 2 + (y?) 2 = 1; this solution corresponds to the vortex points on 
the disk boundary and should be discarded. 

This leaves two families of roots corresponding to 1°: 1 — {x\) 2 + (?/i) 2 = 0, and x\ = 0. 
For family 1° Eq. (j^J) reduces to (y?) 3 — A(y°) 2 + (y®) — A/2 = 0. Solution of this equation 
yields the coordinates x\ and y®, shown in Fig. [2] (left), of stationary locations of the vortex 
points. Note the asymptotic behavior, x\ ~ ±A and y® ~ A as A — > 00. 

It can be seen that there exist stationary configurations of the vortex-antivortex pair, both 
downstream and upstream from the disk; such configurations are symmetric with respect to 
the y-axis. 
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Family ,2° corresponds to the location of the vortex-antivortex pair on the y-axis. Setting 
x = 0, Eq. ® yields (y?) 4 - A(y 1 ) 3 + 2(y°) 2 - 3Ay? + 1 = 0. For y° > 1 (outside the disk) 
the real solution of this equation is shown in Fig. [2] (right). Note the asymptotic behavior 
y® ~ A as A — ► oo. 



IV. APPARENT STABILITY OF STATIONARY POINTS 



The motion of two point vortices whose initial positions are perturbed around the corre- 
sponding stationary values, i.e. 

x J (0)=x J (A) + (Ax,)o, Vj {0) = y%\) + {A Vj ) , (10) 

where j = 1, 2, x°(A) = x^{\), and y%(\) = —yl(X), is governed by the system of four 
Lagrangian equations (jlj) for xx(t), yi(t), x 2 (t), and y 2 ■ Numerical solution of Eqs. (jl]) 
shows that any perturbation leads in general to sweeping of vortex points by the imposed, 
uniform at infinity, flow around the disk, so that Xj(t) — > +oo as t — > oo (the only exception 
being the case of symmetric initial perturbations such that (Axi)o = (^£2)0 an d {Ay 1)0 = 
{Ay2)o, in which case both vortex points move along closed trajectories around stationary 
points). 

Of particular interest would be an analysis of motion of point vortices during the first 
t\ = 70 non-dimensional units of time corresponding to the dimensional duration, t exp = 10 s 
of the experiment (see below Sec. |V|). We will be concerned with the values of A for which 
the vortex points, starting their motion near the stationary points {x*l{\), ±y®{\)), remain 
in the sufficiently close vicinity of their initial positions for at least t± time units. A simple 
a priory estimate for such values of A can be obtained from Eqs. (SI)- (EJ) as follows. At the 
stationary positions, Z\ = z®{\) and z 2 = { z i)* we obviously have fi^{ z i, { Z T)*) = 0- ^ the 
vortex point z\{t) is located initially in the close vicinity of so that |zi(0) — zf \ <^ \z®\, 
the magnitude of its velocity at t — can be estimated, expanding fi around the stationary 
point, as A(A)|zi — z°|, where 

dfl{4: W) 



A{\) 



(11) 



The function A{\) is shown in Fig. [3j 

If, at t — 0, the distance between the vortex Z\ and the stationary point z\ is small, 
the distance traveled by the vortex Z\ during time t\ can be estimated as / ~ A(A)|zi(0) — 



z\\t\. Assuming the distance I to be not larger than few (say, 10) times the initial distance 
between the vortex and the stationary point, we find that, during the first t\ = 70 non- 
dimensional units of time, the vortex point will remain relatively close to its initial position 
provided A(\) < 0.15; this corresponds to A > 3, see Fig. [3J Consequently, on a time 
scale corresponding to t\ non-dimensional units the vortices located sufficiently close to the 
stationary points z°(A) and (z°(A))* will appear, for A > 3, as apparently stable. 

This conclusion is illustrated by the following example. For initial positions defined by 
the perturbations (Axi) = -0.01, (Ayi) = 0.02, (Ax 2 ) = -0.03, and (Ay 2 ) = -0.01, 
the results of numerical calculation of the motion of vortex points are shown in Figs. H] and 
E] (left) for vortices initially located at the rear of the disk, and Figs. [5] and (right) for 
vortices initially positioned in front of the disk. Our calculations show that, for vortices 
initially located both at the rear and in front of the disk, the period of time during which 
the vortex points remain close to their initial locations increases with the non-dimensional 
circulation A. For A = 3 the magnitudes of displacement of vortices from their initial 
locations, Arj = [(xj — x°) 2 + (yj — yj) 2 ] 1 ^ 2 , where j = 1, 2, remain smaller than 0.44 (or 
less than 15% of the distance from the center of cylinder) for t\ = 70 time units. This result 
will be used in Sec. [V] for explanation of the apparent stability of eddies observed in the 
experiment [3j]. (During the same time interval, a more pronounced displacement such that 
Arj ~ 1 occurs already for A = 2.8.) Note that this conclusion is not a consequence of the 
specific choice of initial conditions but remains quantitatively valid for any sufficiently small 
(such that |(A:Cj-)o| < 0.03 and |(Ay 3 -)o| < 0.03) initial perturbation. 

Similar results were also obtained for the motion of point vortices in the vicinity of family 
2 s stationary points located at the y-axis. We found that the point vortices, positioned 
initially near these stationary points, will remain in their close vicinity during the non- 
dimensional time ti = 70 for considerably larger, compared with those for family 1°, values 
of the non-dimensional circulation, i.e. for A > 7; such values of A correspond to the 
stationary points located at y°(A) > 7. Since in the experiment [s| the boundary of the flow 
domain in the ^-direction was at the distance 6.25 cylinder radii, no apparently stationary 
flow structures could be observed corresponding to family U° stationary points. Therefore, 
for the purpose of this work, the further, more detailed, analysis of the motion of point 
vortices in the vicinity of family 2* stationary points would be irrelevant. 

The numerical results described in this Section can explain, entirely from the classical 
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fluid dynamics viewpoint without invoking the mechanism of interaction between the normal 
fluid and superfluid vortices, the apparent stability of the normal eddies observed by Zhang 
and Van Sciver 3| both at the rear and in front of the cylinder in the thermal counterflow, 
see Sec. [V] below. 



V. INTERPRETATION OF EXPERIMENTAL OBSERVATIONS. DISCUSSION 
AND CONCLUSIONS 



In the experiment [3( the velocity field was recorded of a large number of micron-size 
particles injected in turbulent He II thermal counterflow. In classical viscous fluids, provided 
a particle is sufficiently small, the viscous drag force exerted on the particle dominates all 
other forces (such as the inertial and added mass force, Saffman and Magnus lift forces, 
the Basset memory force, etc.), so that the velocity field of a large, dilute ensemble of 
solid particles faithfully maps the fluid flow field. In turbulent He II, due to interactions 
between particles and quantized vortices, particles generally map neither the normal fluid 



nor superfluid, as was demonstrated by recent experimental 



ma 



and theoretical 



studies. However, in Zhang and Van Sciver counterflow experiment |3[ a tangle of 
superfluid vortices was rather dilute, with the mean intervortex distance t ~ 6 /im exceeding 
substantially the particle diameter d p = 1.7 /im. This enables us to assume that close 
encounters between particles and quantized vortices (and, in particular, the events of particle 
trapping on quantized vortex cores) were relatively rare and, therefore, the particulate flow 
field recorded in the experiment [3j maps, at least qualitatively, the velocity field of the 
normal fluid. 

Of particular interest are the apparently stationary normal fluid eddies observed by Zhang 
and Van Sciver both behind and in front of the cylinder. Unlike the familiar Von Karman 
vortex street shed by a cylinder, such structures were never observed in the classical viscous 
flow. In Ref. ^] Zhang and Van Sciver attributed the existence of apparently stationary 
normal eddies to the mutual friction interaction between quantized vortices and the normal 
fluid. 

In this work we argue that, based on the idealized flow model considered above in Sees. HT1 - 
IIVI the experimental results [2I can be interpreted without invoking the mechanism of in- 
teraction between the normal fluid and quantized vortices. It must be emphasized that 



S 



the flow analyzed in Sees. HTllIVI is that of the inviscid fluid, while the normal flow in the 
experiment Q is obviously viscous, so that the following arguments aud estimate, should 
be regarded as qualitative. However, in the experiment [3| the Reynolds numbers defined 
by the diameter of the cylinder were at least of the order of 2 x 10 4 , so that the considered 
inviscid, potential flow can be used as a reasonable approximation to a distribution of the 
Reynolds averaged velocity of the turbulent normal fluid around the cylinder. 

In the experiment described in Ref. the large-scale normal eddies, both downstream 
and upstream of the cylinder, were observed at a distance about 3 cylinder radii at the 
angles ±45° and ±135° to the axis through the center of cylinder in the direction of the 
undisturbed normal flow. For the normal fluid velocity v n « 2.2cm/s, these circulation 
patterns appear as stable for the duration of the experiment t cxp ~ 10 s. 

Our calculation of motion of point vortices in the imposed potential flow around the 
disk showed that there exist stationary locations of point vortices, both downstream and 
upstream of the disk. These locations are unstable: any perturbation of the initial stationary 
positions of point vortices leads, eventually, to sweeping of point vortices away from their 
initial locations. The distance between the stationary points and the center of the disk 
increases with the non-dimensional circulation A = — T/(A7cUa). Also increases with A the 
period of time during which the slightly perturbed vortex points remain in the close vicinity 
of the corresponding stationary points. For A > 3, in the case where the magnitudes of 
initial perturbations of coordinates of stationary positions are smaller than 0.03, we found 
that the vortices will remain in the close vicinity of the stationary points for at least t\ = 70 
non-dimensional units of time. In the case where the undisturbed fluid velocity, U and 
the disk radius, a are identified, respectively, with the normal fluid velocity, v n ~ 22mm/s 
and the cylinder radius, 3.175 mm in the experiment jsj], the corresponding dimensional 
time is t = t\a/U ~ 10 s ~ t exp . Therefore, on a time scale corresponding to the duration 
of experiment jsl the vortices will remain close to the positions of the stationary points 
and, hence, appear as stable. Moreover, for A = 3 the non-dimensional coordinates of the 
stationary points, x° ~ ±3.005 and %f- ~ ±2.834 correspond to the locations, in the cited 
experiment, of the observed apparently stable normal eddies at the angles ±45° and ±135° 
to the axis through the center of the cylinder in the direction of the undisturbed normal 
flow. 

In summary, having considered the motion of the vortex-antivortex pair in the two- 
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dimensional Euler flow around the disk, we found stationary solutions which enabled us 
to interpret the existence of apparently stale normal eddies in the experiment [3J] without 
invoking interactions between the normal fluid and quantized vortices. Nevertheless, the 
mutual friction between quantized vortices and the normal fluid, together with the polar- 
ization of the vortex tangle, might be responsible for the emergence of normal eddies in the 



first p 
Refs. 



ace: the polarized cluster of superfluid vortices would rotate the normal fluid (see e.g. 
) which, in turn, would drag along the tracer particles in the PIV experiment. 
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(On a related problem, it should be mentioned that the results of Hanninen et al. 12j are 
hinting at the possibility of formation of a classical wake of quantized vortices behind an 
oscillating sphere. It is not yet known whether a similar phenomenon occurs in the flow 
past a cylinder, although the experimental results [3j seem to point in this direction.) As- 
suming that, at least within the domains occupied by the circulation cells observed in the 
experiment 3j, the normal fluid and the superfluid are fully interlocked through the mutual 
friction, it is easy to estimate that, for the parameters typical of the cited experiment, the 
value of the non-dimensional circulation, A = 3 corresponds to iV w 2.6 x 10 4 quanta of cir- 
culation, k = 10~ 3 cm 2 /s. Considering a cluster of iV polarized vortices with a cross-section 
of radius a c \, we find that its (polarized) vortex line density would be V w 8.3 x 10 4 cm -2 , 
which is a small, 8% polarization of the total (random) vortex line density, Lq m 10 6 cm -2 



lfl 



111 ] , even a small polariza- 



typical of the experiment [3j . Therefore, as envisaged in Refs. 
tion of the vortex tangle would be sufficient to generate normal circulation patterns which, 
if located near the stationary points of the considered Euler flow, can exist as apparently 
stable for the duration of experiment. 
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FIG. 1: (Color online). Point vortices in the potential flow around the disk. Left: vortex-antivortex 
pair downstream (at the rear) of the disk. Right: vortex-antivortex pair upstream (in front) of the 
disk. 
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FIG. 2: (Color online). Left: coordinates x\ - dashed (red) lines, and y® - solid (black) line of the 
family 1° stationary points in the upper half-plane as functions of the non-dimensional circulation 
A. Right: coordinates of the family stationary vortex point in the upper half-plane (x^ = 0). 
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FIG. 3: (Color online). Function -A(A). The dashed horizontal line corresponds to A = 0.15. 
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FIG. 4: (Color online). Time-dependent coordinates x\{t) (left) and j/i(t) (right) of the vortex 
point initially located in the vicinity of the stationary point (x\{X), y?(A)) in the upper half-plane 
downstream of the disk. Dot-dashed (red) lines: A = 1, x® w 1.191, ps 0.648; dotted (purple): 
A = 1.5, x\ ps 1.554, ps 1.189; dashed (blue): A = 1.8, x° x w 1.828, j/J ~ 1-531; solid (black): 
A = 3, x® ps 3.005, Ui ~ 2.834. Initial positions of the vortex-antivortex pair are defined by 
perturbations (Axi) = -0.01, (Ayi) = 0.02, (Ax 2 )o = -0.03, and (Ay 2 )o = -0.01. 
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FIG. 5: (Color online). Coordinates of the vortex point initially located in front of the disk. Values 
of A, 2/i(A), |x5(A)|, (Axj)o, and (Ayj)o are the same as for Fig. [5j 
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FIG. 6: (Color online). Trajectories of the vortex point. Left: trajectories corresponding to Fig. [H 
Right: first four trajectories, from left to right, correspond to Fig. [SJ the last trajectory corresponds 
to A = 0.5 (x5(0.5) w —1.035, y°(0.5) « 0.266). Inserts show the trajectories for A = 3 during 
t\ = 70 non-dimensional time units. 
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